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We show that complex VT -symmetric photonic lattices can lead to a new class of self-imaging 
Talbot effects. For this to occur, we find that the input field pattern, has to respect specific 
periodicities which are dictated by the symmetries of the system. While at the spontaneous VT- 
symmetry breaking point, the image revivals occur at Talbot lengths governed by the characteristics 
of the passive lattice, at the exact phase it depends on the gain and loss parameter thus allowing 
one to control the imaging process. 



PACS numbers: 42.82.Et, 42.25.Bs, ll.30.Er 

Introduction-The Talbot effect [TJ [2J, a near field 
diffraction phenomenon in which self-imaging of a pe- 
riodic structure illuminated by a quasi-monochromatic 
coherent light periodically replicates at certain imaging 
planes, is an important phenomenon in optics. These 
imaging planes are located at even integer multiples of 
the so-called Talbot distance zt = 2a 2 /A, where a rep- 
resents the spatial period of the pattern and A the light 
wavelength. The simplicity and beauty of Talbot self- 
imaging have attracted the interest of many researchers. 
Such effects find nowadays applications in fields ranging 
from imaging processing and synthesis, photolithography 
[3], and optical testing and metrology [3] to spectrome- 
try and optical computing |5] as well as in electron optics 
and microscopy 6 . Similar processes are encountered 
in other areas of physics involving nonclassical light [7], 
atom optics [3j |8] , Bose-Einstein condensates [9] , coupled 
lasers [10] and waveguide arrays [11]. However all these 
achievements are limited in studying the properties of 
the input beams and using real gratings for imaging. By- 
passing these limitations will not only enrich the conven- 
tional self-imaging research, but also offer new methods 
for imaging technologies. It is therefore extremely desir- 
able to investigate and propose self-imaging architectures 
which incorporate gain or/and loss mechanisms. 

In the present paper we study the Talbot revivals in a 
new setting, namely, a class of active lattices with anti- 
linear symmetries. These structures deliberately exploits 
notions of (generalized) parity (V) and time (T) symme- 
try [T2J H3] in order to achieve new classes of synthetic 
meta-materials that can give rise to altogether new phys- 
ical behavior and novel functionality [141 - H6] . Some of 
these results have been already confirmed and demon- 
strated in a series of recent experimental papers [TSrUT] . 
In classical optics, PT-symmetries can be naturally in- 
corporated [14] via a judicious design that involves the 
combination of delicately balanced amplification and ab- 
sorption regions together with the modulation of the in- 
dex of refraction. In optics, VT— symmetry demands 
that the complex refractive index obeys the condition 
n(r) = n*(— r). It can be shown that these structures 



have a real propagation constant (eigenenergies of the 
paraxial effective Hamiltonian) for some range (the so- 
called exact phase) of the gain and loss coefficient. For 
larger values of this coefficient the system undergoes a 
spontaneous symmetry breaking, corresponding to a tran- 
sition from real to complex spectra (the so-called broken 
phase). The phase transition point, shows all the char- 
acteristics of an exceptional point (EP) singularity. VT- 
synthetic matter can exhibit several intriguing features 
[T4H3T] . These include among others, power oscillations 
and non-reciprocity of light propagation [TH (TS] Q15] , non- 
reciprocal Bloch oscillations [2U] , unidirectional invisibil- 
ity [2S] and a new class of conical diffraction [5T] . In 
the nonlinear domain, such non -reciprocal effects can 
be used to realize a new generation of optical on-chip 
isolators and circulators [22]. Other results include the 
realization of coherent perfect laser-absorber [23l [29] and 
nonlinear switching structures |24j . 

Here, we define conditions which guarantee the exis- 
tence of Talbot self-imaging for a class of active VT- 
symmetric lattices. We find that the non-orthogonality of 
the Floquet-Bloch modes imposed by the non-Hermitian 
nature of the dynamics together with the discreteness of 
the lattice structures imposes strong constraints for the 
appearance of Talbot recurrences. We show that while at 
the spontaneous 'PT-symmetric point the Talbot length 
zt is characterized by the structural characteristics of 
the lattice, in the exact PT-symmetric phase it is con- 
trolled by the gain and loss parameter 7. This allow us 
to have reconfigurable Talbot lengths for the same ini- 
tial pattern. Finally, we discuss possible experimental 
realizations where our predictions can be observed. 

Model- We consider a one-dimensional (ID) array of 
coupled optical waveguides. Each of the waveguides can 
support only one mode, while light is transferred from 
waveguide to waveguide through optical tunneling. The 
array consist of two types of waveguides: type (A) involv- 
ing a gain material whereas type (B) exhibits an equal 
amount of loss. Their arrangement in space is such that 
they form N coupled (A-B) dimers with intra and inter- 
dimer couplings k and c respectively, such that both cou- 
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FIG. 1: (Color online) (a) Photonic lattice structure with 
intra-dimer coupling k and inter-dimer coupling c. Sublattice 
(gain waveguide) A is shown by the red rectangular cuboid 
while sublattice (lossy waveguide) B is shown by green rect- 
angular cuboid. Each dimer is distinguished by the index n. 
(b) Dispersion relations for various 7-values. At 7 = 7pt the 
gap between the two bands dissappear and an exceptional 
point singularity is created. 



plings are of similar (but not the same) size i.e. k ~ c (see 
for example Fig. [2] where k = 1.05c). In the tight binding 
description [35], the diffraction dynamics of the electric 
field amplitude \& n = (a n , b n ) T at the n-th dimer evolves 
according to the following Schrodinger-like equation 



- da n (z) 

; dbj.z) 
1 dz 



ea n (z) + kb n (z) + c& n _i(z) 
e*b n (z) + ka n (z) + ca n+1 (z) 



(1) 



where e = eo + *7 is related to the complex refractive 
index [T3]. Without any loss of generality, we will as- 
sume below that eo = 0, 7 > and c < k [19]. The ef- 
fective Hamiltonian that describes the system commutes 
with an anti-linear operator (in [19] we coined this VdT- 
symmetry) which is related with the local 'PT-symnretry 
of each individual dimer. 

At this point it is beneficial to adopt a momentum rep- 
resentation a n (z) — J dqa q (z) exp(ing) (and simi- 
larly for b n ) where the integral is taken over the Brillouin 
zone —n<q<TT. Because of the translational invariancc 
of the system ([!]), the equations of motion in the Fourier 
representation break up into 2x2 blocks, one for each 
value of momentum q: 



dz 



a q {z) 
b q (z) 



d q (z) 

hi*) 



(2) 



with v q = k + c • e~ lq . The two component wave func- 
tions for different g-values are decoupled thus allowing 
for a simple theoretical description of the system. This 
allows us to perform the evolution in Fourier space and 
then evaluate the spatial representation by a backward 
transformation i.e. 



*„(*) 
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^q(z)e in <*dq. 



(3) 



the n-th dimer in the spatial representation and ip q (z) = 
(a q (z),b q (z)) T is the corresponding Fourier component. 

Dynamics - Substituting in Eq. |2]) the stationary 
form (a n , b n ) T = exp(— i£z)(A, B) T , and requesting non- 
trivial solutions of the resulting stationary problem, i.e., 
(A, B) ^ 0, we obtain the band structure of this diatomic 
V7 system [15] : 



£ ± =±^{k- 



4fcccos 2 (q/2) — 7 2 



(4) 



For 7 = we have two bands of width 2c, centered at 
£ = ±k. In this case, the two bands are separated by 
a gap S = 2(k — c) and the exact VT phase extends 
over a large 7 regime. It follows from Eq. Q that 
when 7 > 7737- = (5/2, the gap disappears and the two 
(real) levels at the "inner" band-edges of the two differ- 
ent bands (corresponding to q — ±tt) become degener- 
ate. The corresponding eigenvectors are also degenerate, 
resulting in an exceptional point (EP) singularity. For 
7 > ypf the spectrum becomes partially complex [19] - 
Below we focus our analysis on the domain 7 < 7737-. 

The eigenvectors associated with the Hamiltonian Eq. 
([2]) are bi-orthogonal, and therefore do not respect the 
standard (Euclidian) orthonormalization condition. As a 
result the conservation of total field intensity is violated 
for any 7 ± 0. Denoting by \R±(q)) = ^(1, £±Mz^ 
the right eigenvectors corresponding to the eigenvalue 
£±(s)j we have that the q— th momentum components 
of any initial excitation can be written as ^(O) = 
Sz=± c i\Ri(l))- The evolved g-field component is 



ipq(z) = ^2c;< 



\Ri(q)) 



(5) 



i=± 



where cf — (Li{q)\tp q (0)) is the expansion coefficient and 
{L[(q)\ is the left eigenvector associated with eigenvalue 
£i(q). The above expansion valid as long as the Hamilto- 
nian H q in Eq. ^ does not have a defective eigenvalue. 
The latter appears at the spontaneous 'FT-symmetric 
point 7737- = k — c (EP) for q — ±7r. The corresponding 
evolved g-field component is then written as: 



4> q =±Tr(z) = (ci + c 2 z)(l, -i) T + c 2 (-i/7, 0) J 



(6) 



where ^> n (z) = (a n (z) , b n (z)) T is the field amplitude for 



Direct substitution of Eqs. ( 5p | into Ecj. Q provides 
the evolution of the field in this system. A note of cau- 
tion is here in order. For the existence of Talbot revivals, 
a necessary condition is that the initial preparation must 
not excite the q — ±7r defective mode. In the opposite 
case, the field increases linearly with the propagation dis- 
tance z (see Eq. ([6])), thus destroying the possibility of 
revivals of any initial pattern. 

Talbot self-imaging - We are now ready to analyze the 
Talbot self-imaging recurrences in the case of the VT- 
symmetric structure of Fig. [T] We recall that in order the 
Talbot effect to occur, the input field distribution should 
be periodic [TT], and thus in general V I'„(0) = ^ n+ ^(0) 
where represents the spatial period of the input field. 
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FIG. 2: (Color online) Talbot intensity "carpets" for period- 
TV input patterns. Length is measured in units of inter-dimer 
coupling c — 1. The intra dimer coupling is k = 1.05. 
(a) periodicity N = 1 with the binary input {1,0,1,0,...} 
and 7 = 0.05. (b) periodicity N — 3 with the binary input 
{1, 1, 1,0,0, 1, 1, 1, 1,0,0, 1, ...} and 7 = ^ VT = 0.05. 



Because of this periodic boundary condition, q can take 
values only from the discrete set 



2m7r 



TV 



-, m = 0, 1,2, ...,JV-1. 



(7) 



Substituting the above constrain in Eq. ^ we get the fol- 
lowing expression for the evolved field at the n-th dimer 



N-l 



* ( nXz)= E cf-e^'^lR^)} (8) 



l=±-m=l 



It is therefore clear that field revivals are possible at inter- 
vals z if £{q m )zT = 2ttv where v is an integer. Therefore 
the ratio of any two eigenvalues £ m = £{q m ) has to be a 
rational number, i.e. 



v /(fc-c) 2 +4fcccos 2 (^)- 7 ; 
^/(fc-c) 2 +4fcccos 2 (^)-7 



a 



(9) 



where a and /3 are relatively prime integers. At the same 
time, revivals in the field intensity are ensured provided 
that (£ m — £^)/ '{£ m i — £^) — a/ (3 where the indices be- 
long to the set {0,1,..., TV — 1} and are taken at least 
three at a time. It is straightforward to show that this 
condition is trivially satisfied for the same set of TV-input 
pattern periodicities as for the fields. 

Next we consider the field Talbot revivals of input pat- 
terns with period TV, at the spontaneous PT-symmetric 
point. To this end, we observe that a direct substi- 
tution of 7 = 7737- in Eq. ([9| for the ratio £ m /£o 
leads to the simple condition cos(mir/N) — a/(3. The 
latter is re-written in terms of the Chebyshev polyno- 
mials which are defined as cos(ma;) = T m (cos(a;)) — 



Y^J2}o c j ( cos ( x )) m 2j ') wn ere [m] represents the integer 
part of m. The Chebyshev coefficients are integer 
numbers and, of importance to our discussion, is the fact 
that the first one is given by = 2 



Given that 



c- are integers, then cos(™) is rational if and only 
if cos(-j^) is rational [TT]. Using the Chebyshev identity 
with m = TV (assuming TV is an odd number), we obtain 
the following polynomial in cos(-^): 



2"- 1 (cob(£))" + 



"[N/2] 



cos(-) + l = (10) 



where we have used the fact that Tjy (cos fj ) — 
cos(Ntt/N) = —1. By applying the rational root the- 
orem one can show that the roots of this polynomial in 
cos(7r/7V) are rational only if N = 1,3. Similar tech- 
niques leads to the fact that for even values of TV the 
only possibility is TV = 2 [11 . However, input pat- 
terns with TV = 2 periodicity, excite the q = ±7r Fourier 
mode, and therefore based on our previous discussion (see 
Eq. (|6| above), have to be excluded. Therefore, strictly 
speaking, discrete Talbot revivals at the spontaneous 
PT-symmetric point are possible only for a finite set of 
periodicities TV — 1,3, where for example, the TV = 1 
case can represent initial patterns {1, 0, 1, • • • ,1,0} or 
{0, 1, 0, 1 • • • , 0, 1} or the more trivial case of a plane wave 
with {1,1, 1, 1 • • • ,1,1}. Some representative intensity 
revivals for TV — 1 and 3 periods are depicted in Fig. [2] 
The Talbot revivals can appear also in the exact phase 
7 < 7737-. Simple inspection of Eq. (9| indicates that 
an initial periodic pattern with periodicity N = 1 (re- 
sulting to eigenvalue index m — in Eq. ([T])) leads 
to a rational value a//3 — 1. In this case the Tal- 
bot length zt depends on the gain and loss parame- 
ter as zt = 27r/£ = 2n/ \Z^j- + ^kc — 7 2 and there- 
fore it varies by changing 7. Such reconfigurable be- 
havior of the Talbot length is characteristic of the ex- 
act phase 7 < tpx and can be found also for the 
TV = 2,3-period input patterns. For TV = 2 (corre- 
sponding to eigenvalue indices m = 0, 1 in Eq. ufh) one 
can show that for fixed k, c and j-pt — k — c such that 
77?r/( fc + c) > Eq. ( [g]) is satisfied provided that 

7 = yj^jyj- — 4fcca 2 / (/3 2 — a 2 ) (we assume that a < 



Similarly for TV 
provided that 7 



3, the Talbot revivals are possible 



kc[l - 4(a/ / 9) 2 ]/[l - (a//?) 2 ] 
where 0.5 < a/ (3 < yl — 3fcc/ (k + c) 2 . In both cases 
the corresponding Talbot length is 7-dependent and it is 
given by the largest period zt — 2n/\£j ~ £i\ ~ 27r/£ 
that results from the eigenvalues involved in the initial 
pattern. Example cases of Talbot self-imaging revivals 
for initial periodic patterns with period TV = 1, 2,3 and 
different 7 values are shown in Figs.[3^-b; Figs.[3j>d; and 
Figs. [3^-f respectively. We see that for the same initial 
preparation, the revivals are controlled by 7 and can oc- 
cur at different Talbot lengths. 

In fact, we can show that larger periods TV > 3 do not 
result in Talbot self-imaging revivals in the exact 'PT- 
symmetric domain. Using Eq.Q, for the |£ m |/|£o| = 
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FIG. 3: (Color online) Talbot intensity "carpets" for period- 
N input field patterns at the exact phase 7 < 7737-. Ev- 
erything is measured in units of inter-dimer coupling c = 1 
while the intra-dimer coupling is k = 4. (a) 7 = 0.1 while 
in (b) 7 = 2. In both cases the input pattern has periodic- 
ity N — 1 and it is chosen to be {1, 0, 1, 0, 1, 0, 1, 0, ■ • • }; (c) 
7 = ^/ll/3 and (d) 7 = \fl. Now the input pattern has 
periodicity N = 2 and it is chosen to be {1, 1, 0, 0, 1, 1, • • • }; 
(e) 7 = 5/^/^7 and (f) 7 = 9/Vl0. The input pattern in 
these cases is {1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, ••• } and has period- 
icity N = 3. Different Talbot lengths zt are observed between 
(a-b) and (c-d) and (e-f). 

a/ (3 and enforcing the constrain that 7 < ^y-pT ~ k — c, 
one obtains the inequality, cos(^) < % which has to be 
satisfied together with the equation Eq. ^ (the equality 
correspond to the case 7 = 7737- discussed above). At 
the same time cos( I ^ L ) has m = 0, • • ■ ,N — 1 roots. By 
applying the intermediate value theorem one finds out 
that this inequality cannot be valid for N > 3. 

Experimental Implementation- We would like finally 
to suggest possible experimental implementations of the 
VT- symmetric waveguide arrays, which will allow for the 
observation of the reconfigurable Talbot effect. The pro- 
posed structures will involve MBE grown quantum wells 
(QW) that will be patterned to form coupled waveguides. 
The basic VT structural element of the array shown in 
Fig. [T] involves two "PT-symmetric sites (dimer) . Such a 
design is desirable because of its simplicity. The dimen- 
sions and index contrast can be such that each waveg- 
uide will be single-moded. For example, for AlGaAs 



structures this can be achieved by a refractive index of 
no = 3.35 operated at 800 nm. Reconfigurable gain can 
be achieved by running an electric current through a Al- 
GaAs/GaAs QW p-n junction. In such structures one 
can easily reach gain and loss values as high as 50cm -1 . 
The two site channels in every dimer will be excited at 
different current levels I\ and I2 so as to establish the 
antisymmetric gain and loss profile that is necessary to 
observe VT optical behavior. In practice this will be 
done provided that current 1\ 3> I2 so as the correspond- 
ing regions underneath see equal amount of gain and loss. 
More specifically, I2 will be relatively small so the associ- 
ated waveguide site will experience material absorption. 
Its sole purpose will be for fine tuning. Given that I\ 
and I2 can be interchanged and adjusted, this will allow 
us to dynamically control the Talbot length zt of these 
PT-symmetric structures. Of course, special considera- 
tion has to be given to the effects of gain and loss on the 
modal index change in these structures (because of the 
Kramers-Kronig relations). 

Finally we comment on the robustness of Talbot re- 
vivals against structural imperfections. For realistic val- 
ues of positional imperfections (up to 5% of the inter- 
dimer coupling) we have confirmed numerically that Tal- 
bot revivals are only slightly distorted. Specifically we 
found that revivals associated with short Talbot lengths 
zt, are essentially unaffected for moderate propagation 
distances z while revivals associated with larger lengths 
Zt are fragile due to the distortion of the delicate balance 
between the mode amplitudes and phases that eventually 
dominate the evolution. 

Conclusions- In conclusion, we have shown that a class 
of PT-symmetric optical lattices, support Talbot self- 
imaging revivals for input patterns with periodicities dic- 
tated by the discreteness of the lattice and the strength 
of gain and loss parameter. Of interest will be to investi- 
gate if Talbot revivals can also occur in higher dimensions 
and in the presence of non-linearity. Our results might 
be applicable to other areas like self-imaging of coupled 
lasers |10j with distributed gain and synchronization of 
PT-symmetric coupled electronic oscillators |16j . 
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